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ABSTRACT 

Measuring the non-Gaussianity of the initial matter density fluctuations may provide powerful 
insights into cosmic inflation and the origin of structures. Although current information on pri- 
mordial non-Gaussianity comes mostly from the Cosmic Microwave Background, obtaining 
constraints from lower-redshift observables will be an important task for future surveys. We 
study the ability of future weak lensing (WL) surveys to constrain primordial non-Gaussianity 
of the local type. We use a large ensemble of simulated WL maps with survey specifications 
relevant to Euclid and LSST. The simulations assume Cold Dark Matter cosmologies that vary 
certain parameters around fiducial values: the non-Gaussianity parameter /nl, the matter den- 
sity parameter Q.^, the amplitude of the matter power spectrum erg, the spectral index of the 
primordial power spectrum Us, and the dark-energy equation-of-state parameter wq. We assess 
the sensitivity of the cosmic shear correlation functions, the third-order aperture mass statis- 
tics, and the abundance of shear peaks to these parameters. We find that each of the considered 
probes provides unmarginalized constraints of A/nl ~ 20 on /nl- Marginalized constraints 
from any individual WL probe are much weaker due to strong correlations between param- 
eters. However, the parameter eiTors can be substantially reduced by combining information 
from different WL probes. Combining all WL probes yields the following marginal (68% 
confidence level) uncertainties: A/nl ~ 50, AQ^ ~ 0.002, Acrg ~ 0.004, Aws ~ 0.007, and 
Awo ~ 0.03. We examine the bias induced by neglecting /nl on the constraints on the other 
parameters. We find erg and wq to be the most affected. Moreover, neglecting non-Gaussianity 
leads to a severe underestimation of the uncertainties in the other cosmological parameters. 
We conclude that a full exploitation of future WL surveys requires a joint analysis of differ- 
ent WL probes. Furthermore, if not taken into account, a non-vanishing level of primordial 
non-Gaussianity will bias the estimated cosmological parameters and uncertainties for future 
surveys. 

Key words: gravitational lensing: weak - cosmological parameters - early Universe - infla- 
tion - large-scale structure - methods: N-body simulations 



1 INTRODUCTION 

According to Big Bang cosmologies, the rich matter structures seen 
today in our Universe have grown from tiny density fluctuations at 
very early times. There are a number of competing theories for the 
physical origins of these primordial density perturbations, e.g. the 
Ekp yrotic model jKhoury et al.||200T| ) or the string gas cosmol- 
ogy ( [Bran denberger et al.'2Q07 ^7The lead ing explanation is the in- 
flationary paradigm (Mukhanov and Chibisov 1981 1. The standard 
single-field slowly-rolling inflationary model provides four funda- 
mental predictions: a nearly flat spatial curvature for spacetime, a 
nearly scale-invariant power spectrum of primordial curvature per- 
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turbations, a nearly Gaussian distribution of primordial density per- 
turbations, and a scale-invariant spectrum of gravitational waves. 
The first three predictions have been tested by observations of the 
cosmic microwave background radiation (CMB) with the Wilkin- 
son Microwave Anisotropy Probe (WMAP, [Komatsu et al.|[2011[ l. 
The discovery of primordial gravitational wave signatures in the 
CMB is one of the main goals for the Planck mission jThe Planck| 
|Collaboration|2006| >. 

Recently much attention has been turned towards using obser- 
vations of the statistical properties of the large-scale structure to 
further falsify models for the origin of density perturbations. For 
example, single field inflation leads to a near-Gaussian distribution 
of fluctuations ^Acquaviva et al.||2003[ [Maldacena 2003 1 . Hence, 
if significant departures from Gaussianity are detected, all single- 
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field inflationary models could be ruled out (I Creminelli and ZaF| 
|darriaga(2004^ . In such models, departures from Gaussianity can 
be characterized by the local expansion of the primordial potential 
perturbations l [Matarrese et al.|2000[ l: 

<b(x) = m + /nl [(f>\x) - {fix))] , (1) 

Here, denotes the Gaussian potential perturbation at space- 
time point X after matter radiation equality, and the parameter /nl 
quantifies the deviation from Gaussianity. In linear theory, typical 
fluctuations are very small with <p ~ 10"'. Thus, for |/nlI < 10^, the 
typical non-Gaussian corrections, are small, too: {fb - (p) I (f) ~ /nl^- 
Currently, the best constraints on local-type primordial non- 
Gaussianity come from the CMB bispectrum in the seven-year 
WMAP data ( |Komatsu et"aL][20TT) : /nl = 32 ± 21 at 68% confi- 
dence level (CL). Though not widely anticipated a few years ago, it 
has recently been shown that probes of structures at lower redshifts 
could also provide competitive constraints, in particular the nonlin- 
ear matter power spectr um ([Se to 2001 , Amara and Refregier 2004 
[Desjacques oi ri^. Tiifepich et al. 2010. Smith et al. 2011), 
the abundance of massive clusters ( .Matarrese et al._2000,,Lo Verde| 
et al.|2008 1, the galaxy bispectrum ([Sefusatti and Komatsu|2007| 



Nishimichi et al.|2010|[Baldauf et al.|2011||Sefusatti et al.|2012| l, 
and the scale-dependence of galaxy clustering l |Dalal et al.|2008( 
[Desjacques eta L 2009 : P illepich et al.|20l"0l l. 

In this paper, we focus on the question: How well can a 
large weak lensing (WL) survey constrain the level of primordial 
non-Gaussianity in the initial conditions? This is of prime con- 
cern given upcoming surveys like the Kilo Degree Survey (KIDS, 
1500 deg^ Kuijken 20 07jl and the Dark Energy Survey (P ES, 



5500 deg , [The Dark Energy Survey Collaboration 



2005 



planned full sky surveys such as those with Euclid (~ 15 000 deg 
[Laureijs et al.[[2011|l or the Large Syn optic Sur vey Telescope 
(LSST,~ 20000dep, [lSST Science Collaborations and LSST 
[Project|2009) . 

There have been various forecasts for the expected error A/nl 
on the non-Gaussianity parameter /nl from future weak lensing 
surveys. For example, Jeong e t aLl ( [2009) predicted that A/nl ~ 50 
could be obtained from galaxy-galaxy lensing on large angular 
scales. [Fedeli and Moscardinl| ( [2010| l found that A/nl < a few tens 
may be possible from a;t-^-analysis of the convergence power spec- 
trum. Using a halo model, Mat uri et al.| ( |201l"| l obtained A/nl ~ 30 
for a combined analysis of the convergence power spectra and the 
abundance of shear peaks in a Euclid-like, survey. Schafer et al.| 
( [2012^ predicted that A/nl ~ 200 could be achieved from ana- 
lyzing convergence bispectra, which could be further reduced to 
A/nl ~ 30 if sufficiently tight priors on all other relevant cosmolog- 
ical parameters are available. Giannantonio et al. (2012) found that 
a combined analysis of weak lensing and galaxy clustering could 
give A/nl < lO. jOguri and Takada[ ( [201 1) estimated A/nl ~ 10 for 
an imaging survey combining cluster counts, the cluster correlation 
function, and stacked WL signals of clusters. 

Whilst this is encouraging, all these studies rely heavily upon 
approximations for the nonlinear matter power spectrum, bispec- 
trum and analytic modelling of shear peak abundances that are not 
well tested or understood in non-Gaussian models. Furthermore, 
the gain in information from combining different probes has not yet 
been fully investigated. Instead, approximations of various degrees 
of plausibility have been made for the cross-covariances between 
different observables. Moreover, most of these forecasts ignore de- 
generacies of /nl with other cosmological parameters. Thus, one 
should be circumspect when interpreting these forecasts. 

In this paper we take a more direct approach to answering 



this question. We use an ensemble of A'-body simulations of struc- 
ture formation in flat Dark-Energy Cold-Dark-Matter (wCDM) cos- 
mologies, through which we perform full gravitational ray-tracing 
to construct a suite of large mock weak lensing maps (e.g. the 
maps for our fiducial cosmological model cover a total area of 
14000deg^). Initial results from these simulations were already 
presented in [Marian et al.[ ( |20I l| l, where we showed that the abun- 
dance of shear peaks was sensitive to /nl, and that one may ob- 
tain an error of A/nl ~ 13 when all other cosmological parame- 
ters are known. In this work we present a much wider study about 
combining information from the cosmic-shear two-point correla- 
tions ^+ and the aperture mass statistics (M^p), and the abun- 
dance of shear peaks (pM- Our approach allows us to fully explore 
the cross-correlations between these different weak lensing probes. 
Furthermore, the initial suite of simulations used in Mari an et al.[ 
l [2011) are augmented with cosmological models that take into ac- 
count variations in other cosmological parameters. Thus, we can 
investigate the joint constraints from WL in the parameter space 
spanned by: the non-Gaussianity parameter /nl, the mean matter 
density flm, normalization erg and spectral index of the matter 
power spectrum, and the equation-of-state (EOS) parameter wq of 
the dark energy (DE). 

Ray-tracing simulations have been used recently by [Pace et al.[ 
( [2011^ and [Fedeli et"aLl (|20TTJ to study the effects of primordial 
non-Gaussianity on basic lensing quantities and the skeleton of 
the convergence field. In contrast to their work, we focus in our 
work on observables more easily accessible in weak lensing sur- 
veys. Moreover, we include realistic levels of noise directly in our 
simulations for the eiTor analysis. Finally, our larger set of simula- 
tions allows us to explore a higher-dimensional parameter space in 
a uniform manner. 

The paper is broken down as follows: In Section|2] we briefly 
discuss the weak lensing cosmological probes that we employ in 
this study. In Section [3] we give an overview of the cosmological 
A'-body simulations and the ray-tracing, and review how we esti- 
mate the observables and their covariance from the simulated data. 
In Section]?] we consider the effects of primordial non-Gaussianity 
on the weak lensing observables. We also discuss the expected sta- 
tistical uncertainties and covariances of the weak lensing observ- 
ables as well as their dependence on the full set of considered cos- 
mological parameters. Furthermore, we present forecasts for un- 
marginalized and marginalized constraints on cosmological param- 
eters, and discuss possible biases arising from neglecting finite non- 
Gaussianity. The paper concludes with a summary and discussion 
in SectionjS] 



2 THEORY 

2.1 Gravitational lensing 

Gravitational lensing is the process whereby photon bundles emit- 
ted from distant light sources are deflected by the gravitational 
field of matter structures along the line of sight. These deflec- 
tions can thus magnify, shift and shear the images of distant galax- 
ies (e.g. [Bartelmann and Schneider|200T| [Schneider et aL|[2006| l. 
The observed image position = (6i,02) of a source at red- 
shift z may thus differ from the source's 'true' angular position 
f) = (J3i{0,z),fi2(0,z))- Image distortions caused by differential de- 
flection can be quantified locally by the distortion matrix 



dp,(0,z) 



'■j=1.2 



"72 + 



~72 ~ ^ 



(2) 
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whose decomposition defines the convergence k{0, z), the complex 
shear y(6,z) = yi(0,z) + iyiiOyZ), and the asymmetry ii)(0,z)- In 
weak lensing, the asymmetry oi is usually small enough to be ne- 
glected (e.g. [Hilbert et al.|2009)|Krause and Hirata|20T0| l. In that 
case, the convergence k and shear y are related by a simple phase 
factor in Fourier space. Thus, in weak lensing, the shear field con- 
tains the same information as the convergence. 

To a good approximation, the convergence is given by a pro- 
jection of the matter density along the line of sight weighted by a 
geometrical factor, which depends on the distances between ob- 
server, deflecting matter structure, and source. The convergence 
field thus bears information on the growth of structures as well as 
on the large-scale geometry of the Universe. 

For a population of sources with known redshift distribution 
p,(z), an effective convergence and shear can be defined by 



k(0) = jdzp,(z)K(0,z) and 
7(0) = Jdzp,(z)y(0,z). 



(3) 
(4) 



Furthermore, the shear y(0) at a position can be decomposed into 
a tangential component yt(0,&) and a cross component y^{0,-&) 
with respect to any given direction & (e.g. [Schneider et al.|2002b) : 

7,(61,1?) = -5l[7(6')e-2"'<*'] and (5) 

7x(e,,?) = -3[7(0)e-^'^""], (6) 

where ip{&) denotes the polar angle of the vector &. 



2.2 Weak lensing observables 

Owing to the fact that we do not know the intrinsic shapes of in- 
dividual background galaxies, the quantities k and y are not direct 
observables. Nevertheless, we can use correlations in the shapes of 
the background galaxies to construct quantities that are observable. 



2.2.1 Cosmic shear correlation functions 

In this study, the first set of lensing observables that we consider 
are the cosmic shear two-point correlation functions (e.g. |Kaiser| 
[T992| >: 



(7) 



where the angle brackets (. . .) denote the expectation for an ensem- 
ble of universes constrained to have a particular set of cosmological 
parameters. The cosmic shear correlation functions f+ and ^_ can 
be directly estimated from the observed ellipticities of the images 
of distant galaxies. 

In weak lensing, the shear correlation functions can be written 
in terms of the convergence power spectrum ^|Kaiser|1992| l: 

r°° dee 

^+/-('?)= —JoMmPM, 

Jq 



(8) 



Jo 2;r 

where J„(x) denotes the Bessel function of the first kind, and P,t(0 
denotes the angular power spectrum of convergence fluctuations. 



2.2.2 Aperture mass statistics 

The aperture mass M^j,(0, &) at position and scale & can be defined 
as a filtered version of convergence or shear field ( [Schneider et al.| 



[T998|: 

Map(6»,!?) = J d^0'U(\0' -0\,-»)k(0') 

= f d-0'Q{\0' -0\,&)y,{0\ff -0), 



(9) 



^ d^0'Q{\0' - 0\, i?)e-^"'<*'-*V(f') 



where 



■ f de'9'U(e', 

Jo 



&) - u(e, &). 



(10) 



Here, we use the third-order aperture mass statistics (M^p)(i?) 
(Mi,p(6', fff) with the filter functions 



1 / eM i_g2 



2(0,1?)= 7Z:^expl 



2&- 



(11a) 
(lib) 



as another lensing probe of cosmology. The statistics {Ml^}(&) can 
also be estimated in a straightforward manner from the observed 
image ellipticities of distant galaxies ( [Schneider et al.|1998^ . 

In the limit of weak lensing, the expectation value of M^^ is 



related to the convergence bispectrum fi^ by (e.g. [Schneider et al.[ 
[2005l >: 



{Ml^)(&) = 2471 1 d{e\ dt'l'l 

Jo J 1 12 Jo 



dif 



(12) 



where f ' = V^^^Tt^^^IZFcos^c, and U(e, &) denotes the (two- 
dimensional) Fourier transform of the filter U(9, i?). 



2.2.3 Shear peak abundance 

The observed shear field can be used to search for massive galaxy 
groups and clusters by applying a filter function matched to their 
expected lensing signal (e.g. Hamana et al. 200' ^. H ere, we con- 
sider the hierarchical filtering method ( M arian et al.[[2009l [201 1| 
[2012| > with aperture-mass filters following the expected lensing 
signal of the Navarro-Frenk-White (NFW) profiles ( [Navarro et al.[ 
[T997) . The observed shear field is convolved with a series of filters 
QmiO', M) of different size, from the largest down to the smallest. 
The size is parametrized by the filter mass M: 



Qje-M)=C{M)y^{0-M\ 



(13) 



where y^Jfl; M) is the shear signal expected for an NFW halo of 
virial mass M at a model redshift z„, = 0.3 (which is roughly where 
halos are most effective lenses for the considered source redshifts). 
These weights (approximately) maximize the signal-to-noise ratio 
for the corresponding peaks ^Marian et al. 2009 2011). The nor- 
malization C(M) is chosen such that the expected filtered signal of 
a halo of mass M at the model redshift Zm equals the halo mass M. 

All maxima in the filtered shear maps are considered as po- 
tential shear peaks. Several selection criteria are used to decide 
whether to regard a maximum indeed as shear peak: i) The max- 
imum must have a signal-to-noise ratio larger than a particular 
threshold, ii) Its amplitude Map must be equal or larger than the 
mass of the filter M used in the smoothing, iii) Its position must 
have a certain minimum separation to other peaks with larger am- 
plitude. 
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A peak will appear with different amplitudes in maps 
smoothed with filters of different size. To assign a unique peak 
amplitude to a peak at position 6, we seek the filter mass M that 
matches the peak amplitude M.^^{0\ M), i.e. we solve the equation 
Map(0;M) = M for M (see [Marian erar][20T2l for details). This 
removes the arbitrariness in the choice of filter scale. Moreover, the 
resulting filter size also maximizes the signal-to-noise ratio 5 IN for 
an NFW halo of virial mass M at the model redshift Zm = 0.3 thanks 
to our particular definition of the filter amplitude. This should also 
(roughly) maximize the sensitivity of the shear peak counts to peaks 
caused by real structures, and thereby to cosmological parameters. 
In the rest of this text, when we refer to the 'mass' of a peak, we 
shall mean the filter mass M for which equation M^p = M is satis- 
fied. 

The number of observed peaks in the filtered fields contains 
information on the abundance of massive structures, and thus on the 
cosmology (e.g. [Marian and Bemstein|2006[ |Dietrich and Hartlap| 
[2010[ l. We define the differential peak abundance ^m(M) at filter 
mass M is by: 



MM) : 



"A51og,o(M)' 



(14) 



where Nnf{M,A) is the number of peaks with peak mass > M in 
the survey area A. This shear peak abundance serves as our third 
cosmological probe. 

In a weak-lensing survey, shear peaks can be readily identified 
and their amplitudes estimated using the observed image elliptici- 
ties of galaxies (see, e.g., Miyazaki et al. 2002 for an early work). 
Unfortunately, to date, there is no satisfactory method to predict 
the abundance of shear peaks for a given cosmology and survey de- 
sign other than numerical simulations of lensing, which we employ 
here. 



2.3 Estimation of cosmological parameters 

We discuss the ability of future weak lensing surveys to constrain 
cosmological parameters within a Bayesian framework. We con- 
sider a collection of model universes jA4(;r)) characterized by a set 
of cosmological parameters n = (n^, . . . ,nf^^y, together with 
a set [d] of A'j -dimensional vectors d = (di, . . . ,dN^y of observa- 
tional data taken in order to constrain the parameters jr. The knowl- 
edge about the probability of having the model universe with pa- 
rameters 7r realized regardless of the outcome of the observations 
is assumed to be quantified by a prior probability p(n). We also 
assume that we know the probability that the observations yield the 
data d, given that the underlying model M has parameters tt. This 
probability is denoted L(7r\d) = p(d\7r), and referred to as the likeli- 
hood. The posterior probability p(7t\d), i.e. the probability of hav- 
ing the parameters jt realized in case the observations have yielded 
the data d, can then be computed by: 



p(n\d) : 



L(n\d)p(n) 



j d'^PTT' L(7r'\d)p(j:')' 



(15) 



2.3.1 The likelihood 

For simplicity, we approximate all relevant probabilities by piece- 
wise constant or multivariate normal distributions. We assume a 
simple multivariate normal likelihood 



-l|l/2 



L(n\d) = 



(27r)"i/ 



- exp 



-i[rf-//(^)]'C,' 



(16) 



where the mean is given by a model prediction fiin), and the covari- 
ance matrix Cd is assumed model-independent. From a large set of 
A'l (Ni » 1) independent realizations rf,- of the data vector for a 
model with parameters n, one can estimate the model prediction fi 
by the sample mean fi, 



1 



//(;r)«=//= 



and the data covariance by the sample covariance Cj 

1 ^' 

Ci " = ^[rf, - fi(n)][d,- - ti(7r)f. 



mated by (Anderson|2003[[Hartlap et aL[2007l l: 



(17) 



(18) 



(Cd)" 



(19) 



1 

which requires N, > N^ + l realizations. 

The maximum-likelihood parameter set n^i{d) for the mea 
sured data d is given by the equation 







dn 



Cd'[rf-//('rmi)]. 



(20) 



A Taylor expansion of the log-likelihood lnL(;r|rf) around the 
maximum-likelihood parameters ;r„,i up to second order in n yields 
the approximation 



L{7:\d) - 



where 



(2;r)~p/2 



exp 



1 



(n - Ttatif C„,, (_7T - n-„,i) 



di^k dpi 

djTj dTTj 



+ (dk-Pk) 



d'p, 



dnjdKj 



(21) 



(22) 



is the inverse of the (maximum-likelihood approximation of the) 
parameter covariance matrix of the likelihood. 

Eq. \2\) reveals that the knowledge to be gained about the 
parameters n from the observations of a particular realization of 
the data d is encoded in the parameter covariance Cn,i(rf) of the 
likelihood. The general information content of the likelihood can 
be quantified by the likelihood Fisher matrix, whose computation 
involves an ensemble average (. . .)j over realizations of the data: 

?^?^r ) ■ (23) 

In our study, where we assume that the likelihood is a multi- 
variate normal distribution with constant parameter covariance and 
model predictions p that are well approximated by affine linear 
functions of the parameters n, the likelihood Fisher matrix can be 
identified with the inverse of the likelihood parameter covariance: 



3n - 3n '''' 
Ignoring all correlations (i.e. approximating Cmi and Cd by 
diagonal matrices), Eq.[24|can be recast into: 



1 



A'd 



hi 



J_Ap^ 



(25) 



where o",,,],, = -yC^Id^ denotes the likelihood standard deviation of 
an element of the parameter vector n, and cTd j = -^Cdj/ denotes 
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the likelihood standard deviation of an element of the data vector 
d. Thus, each term (cr^ldyUj7d;r,) of the sum in Eq. \15\ provides 
a rough measure of how much each datum dj contributes to the 
constraints on the parameter rtj. 

2.3.2 Prior and posterior distributions 

As one case, we consider a flat prior p(7c), i.e. one that is constant 
on some sufficiently large subset 5 c R'^p containing ;r,ni. The pos- 
terior is then a truncated multivariate normal 



and 



p(7T\d) 



(26) 



0, 



where c is a normalization constant, the maximum of the poste- 
rior 7ipo(d) = ;r,ni(rf), and the posterior parameter covariance matrix 

Cpo(rf) = Cml(rf). 

As another case, we consider a Gaussian prior with mean TTp^ 
and covariance Cpj. Then, the posterior is also a multivariate Gaus- 
sian 



p(n\d) 



- exp 



1 



(2;r)"p/2 
with mean 

= [Cp,' + C,;l(rf)]"' [Cp/;rp, + C-J,(d)n^m] , 
and covariance matrix 

Cpo(rf) = [Cp,'+c,;l(rf)f'. 



(27) 



(28) 



(29) 



V^po/ii 



[Cp 



(34) 



From this follows that the posterior conditional variance on a single 
parameter ;r, is given by l/(Cp„')„. 

Furthermore, fixing certain parameters to values differing 
from the unconditional posterior mean may shift the conditional 
posterior mean values of the remaining parameters to values dif- 
ferent from their unconditional posterior mean values. Such shifts 
may be the cause of biases in a parameter estimation. If one or 
more model parameters are wrongly assumed to take values differ- 
ent from their true values, the estimates for the remaining parame- 
ters will also be 'wrong' on average unless there is no correlation 
between the parameters. 



3 METHODS 

We use A'-body simulations of cosmic structure formation along- 
side with ray-tracing to create a large suite of mock weak-lensing 
observations in various CDM cosmologies. From these mock ob- 
servations, we compute the expected signals for the cosmic-shear 
correlation functions, the aperture mass statistics, and the shear 
peak statistics. Furthermore, we estimate the covariance of these 
lensing statistics and their dependence on cosmological parame- 
ters. The results for the parameter dependence and covariance of 
the lensing statistics are then used to forecast parameter errors for 
future large weak-lensing surveys and discuss possible biases aris- 
ing from neglecting non-Gaussianity. 



2.3.3 Conditional and marginal posterior distributions 

For a Gaussian posterior distribution, conditional and marginal dis- 
tributions are also Gaussian. Consider a partition of the parameters 
n = (7ri,JTu) into A'p ; 'interesting' parameters jr-, and A'pu = N^-N^j 
'uninteresting' parameters n^. Similarly, the posterior maximum is 
split into ;rpo = (Tfpo.i, ^po.u)- The posterior covariance matrix and its 
inverse are split into submatrices. 



and C„, 



V^po/11 V po/iu| C^CW 

Here, for example, Cpo.iu denotes the submatrix obtained from the 
covariance matrix Cpo by only keeping the rows of the 'interest- 
ing' parameters and the columns of the 'uninteresting' parameters, 
whereas (Cp^);^, denotes the submatrix obtained by only keeping 
the 'interesting' rows and 'uninteresting' columns of the inverse 
covariance matrix Cp^ . 

The marginal distribution of Tti is then given by 



p(7r,\d) 



(2;r)"p^ 



Ap,,/2 



exp 



^fo.i) Cpo ,i (n, - ;rpo,i) 



(31) 



From this follows that the posterior marginal variance on a single 
parameter rCj is given by Cpo,„ . 

The conditional distribution of is given by 



p(7Ti\d,7Tu) - 



where 



c: 



-l|l/2 



(2;r)"p. 



exp 



1 



n.MnJ = Kp^Ad) - [(Cp^)iJ"'(Cp„'),„ [n,, - ;rpo,„(rf)] 
= ^pa\{d) + Cpo,iu(Cpo_„i,) ' [;r„ - ;rpo,u(rf)] , 



(32) 



(33) 



3.1 N-body simulations 

We use two sets of A'-body simulations of structure formation in 
flat CDM cosmologies for our work. To study the dependence 
of weak lensing observables on primordial non-Gaussianity, we 
use the (non-)Gaussian simulations described in |Desjacques et dl] 
1 ,2009. 1 and ,Marian et al.|p01l"| l. These assume a flat ACDM cos- 
mology whose parameters are compatible with WMAP 5-year data 
l |Komatsu et al.|2009t : h = 0.7, n„ = 0.279, fib = 0.0462, Dde = 
0.721, Wo = -1, a-8 = 0.81, n, = 0.96, and /nl e {-100,0, 100). 
The simulations employ Gadget-2 ("Springel 2005 1 to follow 1024^ 
particles in a cubic volume of 1600A"' Mpc side length with a force 
softening length of 40/;"' kpc. 

There are six (non-)Gaussian simulations for each value of 
/nl. The initial conditions of the simulations are matched between 
the different /nl values, which reduces the influence of cosmic vari- 
ance on the relative values of measured weak-lensing observations 
for difi'erent /nl. We choose the six simulations with /nl = to 
represent our fiducial cosmology. 

To study the dependence of weak lensing observables on a 
number of other parameters, we make use of the zHORIZON sim- 
ulations by |Smith| p009^ . These simulations use Gadget-2 with 
750' particles to follow the structure formation in a cubic region 
of 1500/!"' Mpc side length in flat CDM cosmologies with param- 
eters around the best-fitting WMAP 5-year values. 

There are eight zHORIZON simulations for a flat CDM cos- 
mology with h = 0.7, n,„ = 0.25, fit = 0.04, Dde = 0.75, 
Wo = -1, o"8 = 0.8, iis = 1, and /nl = 0. Furthermore, there 
are 32 simulations that vary one of the above parameters, four 
each for £1„ 6 {0.2,0.3), erg 6 {0.7,0.9), n, e {0.95,1.05), and 
Wo 6 {0.8, 1.2). Together with the (non-)Gaussian simulations. 
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these allow us to study the change of observables with varying /nl. 
flm, tTs^ 'Js, or Wq. 



3.2 Lensing simulations 

For the lensing simulations, we closely follow the method of 
|Hilbert et al.| ( |2009| ), which employs a Multiple-Lens-Plane ray- 
tracing algorithm (e.g. [Blandford and Nar ayan 1986) |Jain et aT] 
|2000| l to calculate the light deflection by matter structures in the 
simulations. The past lightcone of a fiducial observer in a simula- 
tion is partitioned into redshift slices with about 120/;"' Mpc co- 
moving line-of-sight thickness (assuming the flat-sky approxima- 
tion). Each slice is then filled with the matter structures of the sim- 
ulation snapshot closest in redshift to the mean redshift of the slice, 
and the slice content is projected onto a lens plane. Light rays are 
then traced back from the observer through these lens planes, at 
which the light rays are deflected. The deflection angles (and their 
derivatives) at the lens planes are calculated from the projected mat- 
ter distribution by a mesh method with a resolution of 20/i"' kpc 
comoving. 

For each simulation, we generate 16 (quasi-)independent 
fields of view by choosing different observer positions within the 
simulation box. Each field has an area of 12 X 12deg^ and IS cov- 
ered by a regular mesh of 4096^ pixels in the image plane, which 
provides an angular resolution of 10 arcsec. Thus we obtain 96 sim- 
ulated survey fields covering a total area of 14000deg^ for each 
/NL-value from the (non-)Gaussian simulations. Furthermore, we 
have 128 fields and a total area of 18 000deg" for the fiducial cos- 
mology of the zHORIZON simulations, and 9 000 deg" for the other 
cosmologies of the zHORIZON simulations. 

For each pixel in the simulated fields, a light ray is traced back 
through the planes up to source redshift z = 4. The ray distortions 
along the rays are then used to calculate the effective convergence 
and shear for a population of sources with median redshift Zmedian = 
0.9 and redshift distribution jBrainerd et al.|1996t 



Pziz) 



3z^ 

— exp 



2z, 



Zo 



I 
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O.TlZmc. 



(35) 



This redshift distribution is similar to those expected for the 
plaime d lensing surveys with Euclid ([Lau reijs et al.] 2011^ o r the 
LSST jLSST Science Collaborations andLSST Project|2009^ . 

The resulting simulated shear and convergence maps are free 
of noise due to the finite number of source galaxies and their in- 
trinsic ellipticities. For a realistic description of the measurement 
uncertainties, we also need to take this shape noise into account. 
We thus create shear maps with shape noise by adding Gaussian 
noise with standard deviation cTpi,^ = cr^. / ^JruA^ to each pixel 
(with pixel area Apix ~ 10^ arcsec^) of the simulated shear fields. 
We assume an intrinsic galaxy ellipticity distribution with standard 
deviation a^. =0.3 for each ellipticity component, and a source 
galaxy density = 40arcmin"^. 



3.3 Estimation of the model predictions 

We use the simulated lensing fields to compute the expected signals 
for the cosmic-shear correlation functions f+()?) in bins of separa- 
tion 1?, the third-order aperture mass statistics {M^^iff) as a func- 
tion of filter scale &, and the shear peak abundance (pM in bins of 
filter mass M. 

For each field, we estimate the shear correlation functions 
from the shape-noise free shear maps employing the Fast Fourier 



Transform (EFT) technique described in |Hilbert et al.|p01 \ ). The 
two-dimensional auto- and cross-correlations of the shear compo- 
nents 7i and 72 are computed employing EFTs. The resulting two- 
dimensional correlations (7,y/)(i?) are then multiplied by direction- 
dependent phase factors depending on the direction of & and binned 
according to separation to obtain the binned cosmic shear corre- 
lations. 

Estimates for the aperture mass statistics are computed from 
the shape-noise free convergence maps of the fields. The conver- 
gence maps are convolved with the filter \\\&\ of a given scale 1? 
by using EFTs and exploiting the convolution theorem. The pre- 
diction for the aperture statistics {Ml^){-ff) is then computed by a 
simple average of the cubed values of the filtered convergence field 
M^p(0, 0), whereby points closer than 4& to a field boundary are 
discarded. The procedure is carried out for the different filter scales 

Estimates for the shear peak abundance are obtained from the 
noisy shear maps using the hierarchical algorithm described in Sec- 
tion 2.2.3 ' We employ a series of 12 filters logarithmically span- 
ning the range M = 8.85 x lO'^/i"' Mq to 2 x lO'^/r' Mq, with a 
signal-to-noise detection threshold (5 /N)^i„ = 2.6. Linear interpo- 
lation is used to estimate the filter mass M with the best matching 
amplitude for each peak. Finally, peaks too close to larger peaks 
(i.e. within the projected virial radius of an NEW halo at the model 
redshift with the same amplitude as the larger peak) are discarded, 
and the remaining peaks are binned according to their peak mass 
M. Note that we do not use any information about any actual halos 
that create some of the peaks. 

For each cosmology, we calculate model predictions for the 
considered lensing statistics by a simple average of the statistics 
estimated in each of the 64 to 128 simulated fields of that cosmol- 
ogy according to Eq. ( |17^ . The model predictions for the differ- 
ent lensing statistics and scales are then combined into one vector 
fi = iUthlU&i), ■ ■ . ,^-('?i), . . .,(M3p)(^i), . . ■,'Pm(M,\ . . .). 

We assume for simplicity that the model predictions // 
are linear functions of the cosmological parameters n = 
(On,, (Tg, h'o,/nl) in the range covered by our simulations. We 
compute the derivatives of the model predictions w.r.t. the cosmo- 
logical parameters by a central difference: 



1 .1 



(36) 



Here, ;r'^-^'' and are two sets of parameters that differ only 

in the y-th component. Our particular choices of parameters for the 
simulations (see Sec. |3.1^ ensures that we have predictions for two 
suitable sets of parameters available for all partial derivatives. We 
also check that using forward or backward differences instead cen- 
tral differences does not yield significant differences in the esti- 
mated derivatives. 



3.4 Estimation of the measurement uncertainties and data 
covariances 

We estimate the covariance of the lensing statistics from the sample 
covariance of the measured statistics in the 96 12 x 12 deg^ simu- 



' Unlike the predictions for the shear correlations or the aperture statistics, 
the model predictions for the peak abundance depend on the level of shape 
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lated fields created from our Gaussian simulations [^To obtain real- 
istic covariance estimates, we measure the shear correlation func- 
tions, the aperture mass statistics, and the shear peak functions from 
the shear maps with shape noise. 

For the noisy measurements of the shear correlations functions 
we use the same method as for the model predictions, but with 
the noisy shear maps as input. The noisy aperture mass maps for 
measuring (M^p) are obtained by filtering the noisy complex shear 
maps with the filter \\Vo) modified by a phase factor to select the 
tangential shear component. The shear peak counts are obtained 
from the same noisy shear maps using the hierarchical filtering. The 
resulting measurements of {M\^, and (f>M in bins of separation 
or scale obtained from each noisy shear field r are arranged into a 
noisy data vector rf,- = . . .)'. 

The data covariance matrices and inverse data covariance ma- 
trices are calculated from the sample covariance of the noisy mea- 
surements in each simulated field by using Eqs. (18) and The 
resulting matrices describe the expected data covariance for one 
field of 144 deg^. We rescale the matrices to describe covariances 
of a survey with 125 independent 144 deg^ fields, i.e. with a to- 
tal area of ISOOOdeg-. With that area, a median source redshift 
^median - 0.9, a source galaxy density = 40arcmin"^, and a stan- 
dard deviation cr^_ = 0.3 per ellipticity component, the resulting 
covariances are representative of a large future 'all-sky' weak lens- 
ing survey similar to those planned for Euclid or the LSST. 

The covariances calculated in that way probably underesti- 
mate the statistical errors for such large surveys. Since our sim- 
ulated fields are created from only six independent simulations 
and therefore are partially correlated, the sample covariance of the 
fields might underestimate the true covariance expected for such 
fields. Furthermore, due to large-scale cosmic variance, the data 
measured in a survey covering half the sky is more correlated than 
data from a hypothetical survey with the same total area but dis- 
tributed in many independent fields. 

We are also interested in the statistical uncertainty of our 
model predictions and the relative importance of shape noise and 
cosmic variance. If only shape noise and cosmic variance are con- 
sidered as sources of noise, the covariance in the measured shear 
correlation functions and aperture statistics can be split into sev- 
eral parts: a pure shape noise part, a pure cosmic variance part, and 
mixed noise parts due to a coupling between shape noise and shear 
correlations (e.g. [Schneider et al.|2002a| |KiIbinger and Schneider| 
|2005j ). We estimate the cosmic variance contribution to the cosmic 
shear correlation functions and the third-order aperture statistics 
from the sample covariance of the shape noise-free simulations. 
The expected pure shape noise and mixed noise contributions are 
computed from the assumed source galaxy density, intrinsic ellip- 
ticity variance, and measured shear correlations. 



4 RESULTS 

In the following, we discuss the effects of primordial non- 
Gaussianity of the local type on weak lensing observables, and 
how these effects can be used in large surveys to constrain the non- 
Gaussianity parameter /nl and other cosmological parameters. In 
particular, we present the separate and joint parameter constraints 

^ As a check, we also compute the covariances from the simulated fields 
of the non-Gaussian simulations, and from the fields of the fiducial model 
of the zHORIZON simulations. The resulting covariances are very similar 
to those estimated from the Gaussian simulations. 



obtained from cosmic shear correlations, third-order aperture statis- 
tics, and shear peak statistics. Furthermore, we discuss the biases 
in other cosmological parameters induced when a non-zero non- 
Gaussianity is neglected in the parameter inference from weak lens- 
ing data. 

As indicated in the outline of our simulations, we discuss re- 
sults expected from a Euclid- or LSST-like lensing survey with an 
area of 18 000 deg^, a source galaxy population with a median red- 
shift Zmcdian = 0.9, a uumbcr density = 40arcmin"^, and a stan- 
dard deviation a^. = 0.3 per intrinsic ellipticity component. We 
assume the correlation functions ^±(&) are measured in 10 logarith- 
mically spaced bins in the range 0.5 arcmin < i? < 120 arcmin. The 
aperture mass statistics {Ml^)(-&} is measured for 10 logarithmically 
spaced filter scales & in the range 0.5 arcmin < j? < 60 arcmin. The 
peak abundance ^^(M) is measured in 10 logarithmically spaced 
bins in the range 10"* < M < 1.3 x 10'^ Mq. 

We only consider uncertainties in the measurements due to 
shape noise (i.e. due to the unknown intrinsic shapes of individual 
galaxies) and cosmic variance (i.e. due to the variation in structures 
between particular portions of the Universe), which set a funda- 
mental limit on the accuracy of the results from a weak lensing sur- 
vey. We do not account for other contribution to the error budget 
(e.g. from uncertainties in the model predictions or biases related 
to imperfect galaxy shape and redshift measurements in observa- 
tions). 

4.1 Effects of primordial non-Gaussianity on lensing 
observables 

First, we study how primordial non-Gaussianity changes the ob- 
served shear correlation functions, aperture statistics, and peak 
statistics compared to the Gaussian case. The changes in the lens- 
ing statistics due to non-Gaussianity are usually small compared to 
the dynamic range in which the signals vary in the considered range 
of angular scales or masses. We thus discuss relative differences in 
the signals for different /nl- 

Figure [T] shows the relative differences in the cosmic shear 
correlation functions f+ and ^_ between the Gaussian case /nl = 
and the non-Gaussian case /nl t 0. The error bars indicate the sta- 
tistical errors due to the intrinsic galaxy image shapes and cosmic 
variance expected for a large future survey (see Section [42| for a 
more detailed discussion). 

On scales & < I deg, the shear correlation function ^+()?) is 
larger for larger /nl. The largest relative difference is visible on 
scales J? ~ I arcmin. The difference is only about ±2% for /nl = 
±100, but this is large enough to be detectable in future surveys, 
since the errors are < 1% at these separations, as shown in the 
next section. The results for ^_ are very similar, but with the largest 
relative differences shifted to larger scales ~ 10 arcmin. 

Both ^+ and ^_ are weighted projections of the three- 
dimensional matter correlation, with ^+(&) probing larger 3-D 
scales than ^_(i?) at the same angular scale &. At redshifts z ~ 0.3, 
where lensing is most efficient, the angular scales where the shear 
correlations show the largest response to /nl correspond to 3-D 
scales of a few hundred kpc. This result is consistent with predic- 
tions from a halo model ( Smith et al.||201l} . In the presence of 
primordial non-Gaussianity, both the abundance and density pro- 
files of massive haloes are slightly increased/decreased for posi- 
tive/negative /nl- The relative contribution of these massive ha- 
los to the matter correlation is largest at the group and cluster 
scales. On smaller length scales, the signal becomes dominated 
by the lower mass halos, for which the effects of primordial non- 
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Figure 1. Cosmic sheai' coirelations (■&) as functions of separation 
1? for different values of non-Gaussianity /nl. Points show the relative dif- 
ference between /nl and /nl = 0. Error bars indicate the statistical 
uncertainties due to shape noise and cosmic variance expected for large fu- 
ture surveys. 



Gaussianity are significantly reduced. On larger scales, two-halo 
terms, which are also less affected by primordial non-Gaussianity, 
dominate the matter correlation. This signature is qualitatively mir- 
rored in Figure [T] 

Figure |2] shows the differences in the aperture statistics 
(M;Jp)()?). Similar to the shear correlation functions, the third- 
order aperture mass statistics is larger for larger primordial non- 
Gaussianity. The relative difference in {Ml^iff) is about ±5% for 
_/nl± 100, which is clearly larger than the statistical uncertainties of 
about 1% expected for scales i? < 5 arcmin in large future surveys, 
as shown in the next section. 

The impact of primordial non-Gaussianity on the shear peak 
abundance (jy^iM) is shown in Fig. [3] The peak abundance in- 
creases with increasing /nl- Similar results were presented in our 
previous study ( [Marian et al.||20fT| . The change is stronger for 
larger peak masses M. For these larger masses, however, also the 
noise becomes larger 



4.2 Error analysis 

As discussed in Section [J!4] the uncertainties in the measured shear 
correlation functions can be split into a pure shape noise part, a pure 
cosmic variance part, and a mixed noise part due to a coupling of 
shape noise and shear correlations. In Fig. |4] these contributions 
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Figure 2. Third-order aperture mass statistics {M\^)^^^^\p) as function of 
filter scale for different values of non-Gaussianity /nl- Points show the 
relative difference between /nl and /nl = 0. Error bars on indicate the 
statistical uncertainties due to shape noise and cosmic variance expected for 
large future surveys. 




Figure 3. Shear peak abundance (p (M) as function of peak mass M for 
different values of non-Gaussianity /nl. Points show the relative difference 
between /nl ^ and /nl = 0. Error bars indicate the statistical uncer- 
tainties due to shape noise and cosmic variance expected for large future 
surveys. 



are shown together with the full statistical uncertainties for our as- 
sumed survey. The cosmic variance dominates the error on ^+ on 
all considered scales (cf., e.g., [Beynon et al.||20l"0} . The cosmic 
variance sets a lower limit of ~ 0.5% on the uncertainty in 
for separations i? ~ 1 arcmin, and of > 1% on i? > 10 arcmin. 
Shape noise is the dominant contribution only for ^_ on scales 
!? < 2 arcmin. On larger scales, the cosmic variance and mixed 
noise dominate the error on 

As previously mentioned, our covariance results based on the 
simulated 12 x 12deg" fields might underestimate the covariance of 
^^i^) in a full-sky survey, in particular for large separations i?. To 
check this, we also compute the expected full-sky covariance em- 
ploying the normal and the log-normal approximation ( [Schneider| 
|et al.|2002b{ [Hilbert et al.|2011| l. The results from these analytic 
approximations indicate that our simulation results do not underes- 
timate the covariance except for f + on large scales i? > 1 deg, where 
the difference is ~ 20%. 
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Figure 4. Expected contributions to the relative errors (symbols, connected 
by lines for visual guidance) on the estimates of the shear correlation func- 
tions ^±(ff) in bins of separation ??: full en'or (squares), intrinsic ellipticity 
noise (diamonds), mixed noise (stars), and cosmic variance (triangles). 




Figure 6. Expected relative errors (symbols, connected by lines for visual 
guidance) on the estimates of the shear peak abundance ipM(M) for various 
peak masses M. 

Figure [5] shows that the cosmic variance is also the dominant 
contribution to the error on the aperture statistic (M;Jp)(i?) on scales 
& > I arcmin. Only for smaller filter scales, there is a significant 
contribution from galaxy shape noise and mixed noise. 

The estimated errors for the shear peak abundance are shown 
in Figure[6] The errors are ~ 1% for the bins with peak masses of 
a few 10'^"' Mq or smaller, but rapidly increase to tens of percent 
for larger peak masses. There is no easy way to split the errors into 
contributions from shape noise or cosmic variance. However, by 
comparing results from fields with the same cosmic structure but 
different shape noise realizations, we find that shape noise makes a 
substantial contribution to the statistical errors on the peak counts 
even for intermediate and large peak masses. A more thorough un- 
derstanding of the noise properties of the shear peak statistics is 
certainly worthwile in future. 
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Figure 5. Expected contributions to the relative errors on the estimates 
of the third-order aperture mass statistics (M;Jp)(i9) for various filter scales 
t?: full error (squares), intrinsic ellipticity noise (diamonds), mixed noise 
(stars), and cosmic variance (triangles). 



4.3 The covariance of the lensing statistics 

The error analysis in the previous section was restricted to the au- 
tovariance of the lensing observables ^±(9), {Ml^}(ff), and (puiM). 
Our lensing simulations allow us to estimate the full covariance be- 
tween these observables, quantified by the covariance matrix Cj. 
For an analysis of the correlation between the observables, we con- 
sider the correlation matrix 



■<i.ii I 



(37) 



The correlation matrix is shown in Fig. |7] All f+-bins with 
similar separations ■& are strongly correlated with each other. The 
correlation among the ^_-bins and among the (M;Jp)-bins is slightly 
smaller. There is also a substantial correlation between and 
(M^p). Such large correlations are expected for surveys like the one 
considered here, where the cosmic variance is the dominant source 
of uncertainties. 

In contrast, there is little correlation between different bins of 
0M- The correlation between (pu and the other lensing observables 
is not very large either. This indicates that by combining the shear 
peak statistics with the other probes, one can obtain substantially 
more information than one could get from either shear peaks or, 
e.g., cosmic shear correlations alone (as advocated by |Dietrich and| 
|Hartlap|2010) . 
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Figure 7. CoiTelation matrix quantifying tlie correlation between the entries of the full data vector, containing the shear correlations functions ^±(6) measured 
in 10 bins of pair separation 9, the aperture statistics (M^p)()?) measured for 10 filter scales 6, and the shear peak abundance (pu(M) measured in 10 bins of 
peak mass M. Within each class of observables. filter masses, scales or separations increase from left to right and from top to bottom. Darker colors indicate 
larger absolute correlation. 



4.4 Derivatives of the lensing observables 

In order to translate the measurement uncertainties and correlations 
in the lensing observables into uncertainties and correlations in the 
cosmological parameters, we need to know how the observables re- 
spond to changes in the parameters. As apparent from Eq. ([24j, the 
derivatives of the model prediction with respect to the parameters 
are sufficient for this. 

The derivatives we estimate from our simulations via finite 
difference are shown in Fig. [8] The Figure also shows the statistical 
uncertainties on our results for the derivatives stemming from our 
limited number of realizations. In most cases, the derivatives bear 
a relatively small statistical uncertainty. Only for {Ml^} and 
the uncertainties in the derivatives reach a magnitude similar to the 
derivatives themselves on larger aperture and mass scales. 

In order to understand how much each observable contributes 
to the parameter constraints, and how much the uncertainties in the 
derivatives might affect the estimated parameter errors, we consider 
the individual contributions to the sum in Eq. |25j. These contribu- 
tions, the derivatives rescaled by the expected statistical errors on 
the observables due to shape noise and cosmic variance, can be con- 
sidered as measures of differential signal-to-noise. They are shown 
inFig.|9] 

The scales with the largest contributions to constraints on /nl 
from ^+ and f_ are ~ 1 arcmin and ~ lOarcmin, which are the 
scales where the relative changes in ^+ and f_ with changing /nl 
are largest. For {Mj^}, the largest contributions are from scales 
~ 1 arcmin, where the relative errors are smallest. For iJm, the 



largest contributions to constraints on /nl are not from large peak 
masses, where the relative change with changing /nl is largest, but 
from intermediate peak masses due to the much smaller noise|^ 

The results for the contributions to the constraints on the other 
cosmological parameters are similar. In general, most constraints 
come from small and intermediate scales, with the largest angular 
and mass scales contributing relatively little. In particular, the peak 
counts in the small mass range appear to provide the largest contri- 
butions to constraints on the DE EOS parameter wq. 



4.5 Cosmological constraints 

In this section, we discuss the cosmological constraints expected 
from measuring and analyzing the cosmic shear correlation func- 
tions the third-order aperture statistics (M^p)()?), and the 
shear peak statistics 0m(M) in a large future WL survey. 

We restrict the analysis to flat CDM cosmologies with a dark 
energy with a time-independent equation of state and local-type 
primordial non-Gaussianity, which is the class of models covered 
by our suite of simulations. The considered parameters are the non- 
Gaussianity parameter /nl, the mean matter density Q.^, the matter 
power spectrum normalization erg and spectral index n^, and the 
dark-energy equation-of-state parameter wq. Here, we concentrate 
on the conditional and marginal parameter constraints, as well as 



^ Note however that this simple consideration ignores correlations between 
bins. 
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Figure 8. Derivatives of tlie cosmic shear correlations the aperture statistics (M^p), and the shear peak abundance (f>M with respect to the non-Gaussianity 
parameter /nl, the mean matter density Cl^, the matter power spectrum normalization erg and the spectral index rig, and the equation-of-state parameter wq of 
the dark energy, as estimated from our simulations (points with error bars indicating uncertainties on the estimate). 



possible biases due to ignoring non- vanishing non-Gaussianity. The 
full posterior parameter covariances for the considered probes are 
listed in Appendix [B] 

4.5.1 Conditional constraints on cosmological parameters 

The conditional posterior uncertainties expected from our model 
survey on individual cosmological parameters, when flat priors are 
assumed, are listed in Table[T] When all other cosmological param- 
eters are known, the cosmic shear correlations the third-order 
aperture statistics (M^p), and the shear peak abundance each 
yield 1-sigma constraints of A/ml ~ 20, with the smallest error 
A/nl ~ 17 from 0m. When all WL probes are combined, the con- 
ditional error A/nl « 14. 

The cosmic shear correlation functions and shear peak abun- 
dance also provide similar conditional constraints on the other 
cosmological parameters f2,„, erg, ;is, and wq. The parameter con- 
straints from the aperture statistics are generally weaker than those 
from ^± or (pM. By combining two or all three of the considered 



Table 1. Conditional constraints on cosmological parameters from a large 
future weak tensing survey: estimated conditional posterior standard devi- 
ations (Tc of the non-Gaussianity parameter /nl, the mean matter density 
flni, the normalization erg and the spectral index of the matter power 
spectrum, and the equation-of-state parameter wq of the dark energy, ob- 
tained when using the cosmic shear coiTelations f + , the aperture statistics 
(M'p), or the shear peak abundance i^m either alone or combined with each 
other. 
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Figure 9. Derivatives of the cosmic sliear correlations f+ (1st and 2nd column), tlie aperture statistics (M^p) (3rd column), and tlie sliear peak abundance i^m 
(4th column) with respect to the non-Gaussianity parameter /nl (top row), the mean matter density (2nd row), the matter power spectrum normalization 
D"8 (3rd row) and the spectral index (4th row), and the equation-of-state parameter vvo of the dark energy (5th row) relative to the statistical eiTors (Tj on the 
observables predicted for a large future survey (points with eiTor bars indicating uncertainties stemming from our noisy estimates of the derivatives). 



WL probes, the errors on fi,^, erg, and wq can also be improved. 
The relative conditional constraints on these parameters expected 
from a combined analysis of all the WL probes are ~ 0.1% - 1%, 
as shown in Table [T] 

The joint conditional constraints for pairs of cosmological pa- 
rameters are illustrated in Fig. [To] For all considered WL probes, 
/nl shows strong (anti-)correlations with the other considered cos- 
mological parameters. (Anti-)correlations are also visible between 



the other parameters, e.g. between erg and £1,^. The constraints from 
combining the WL probes benefit substantially from the comple- 
mentary constraints provided by the different WL probes. Com- 
bining the information from the different probes also reduces the 
parameter correlations. 
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Figure 10. Joint conditional constraints on pairs of cosmological parameters from a large future weak lensing survey. Considered parameters are the non- 
Gaussianity parameter /nl, the mean matter density Hm, the normalization erg and the spectral index «s of the matter power spectrum, and the equation-of-state 
parameter wo of the dark energy. Shown are the estimated 68% confidence contours obtained from the cosmic shear correlations f+ (dashed lines), the aperture 
statistics {M^ ) (dotted lines), the shear peak abundance 0m (dash-dotted lines), all of these WL probes together (solid lines). 



Table 2. Marginal constraints on cosmological parameters from a large fu- 
ture weak lensing survey: estimated posterior standard deviations (Tp„ of the 
non-Gaussianity parameter /ml, the mean matter density n,,,, the normal- 
ization (Tg and the spectral index of the matter power spectrum, and the 
equation-of-state parameter wq of the dark energy, obtained when using the 
cosmic shear correlations the aperture statistics (A/^p), or the shear peak 
abundance ipM either alone, all together, or all together in conjunction with 
results from a CMB experiment like Planck. 
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4.5.2 Marginalized constraints on cosmological parameters 

The marginal posterior standard deviations on the cosmological pa- 
rameters, assuming flat priors, are shown in Table |2] The marginal 
constraints on /nl from (Ml^) alone are so large (A/nl ~ 500) 
that they are practically useless due to the relatively large uncer- 
tainties in the other cosmological parameters and their correlation 
with /ml. The marginal errors on ^nl from f± (A/nl ~ 200) or 
(pM (A/nl ~ 70) alone are much smaller than the constraints from 
(M^p) thanks to the thighter constraints on the other cosmologi- 
cal parameters. All WL probes combined yield a marginal error 
A/nl ~ 50, which is somewhat larger than current constraints from 
CMB bispectra (A/nl ~ 20). 

The marginal errors from either the shear correlations, the 
aperture statistics, or the peak counts on the other cosmological 
parameters are ~ 1% - 20% of the fiducial values. However, com- 
bining the information from cosmic shear and the shear peak abun- 
dance yields much more competitive marginal constraints, with un- 
certainties much smaller than those from any of the individual WL 
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Figure 11. Joint marginal constraints on pairs of cosmological parameters from a large future weak lensing survey. Considered parameters are the non- 
Gaussianity parameter /nl, the mean matter density Hm, the normalization erg and the spectral index of the matter power spectrum, and the equation-of-state 
parameter wo of the dark energy. Shown are the estimated 68% confidence contours obtained from the cosmic shear correlations f+ (dashed lines), the aperture 
statistics (M^p) (dotted lines), the shear peak abundance <I>m (dash-dotted lines), all of these WL probes together (solid lines), or all WL probes combined with 
results from a Planck-\\ks CMB experiment. 



probes. This is in accordance with findings by [Dietrich and Hart-| 
|lap| ( |2010l l, whose study however was limited to n„ and erg. An 
additional improvement can be obtained by also adding informa- 
tion from the aperture statistics. The great improvement reached by 
combining the three WL probes compared to only analyzing any 
single WL probe is also evident in the joint marginal constraints 
shown Fig.|l 1| 

Table [2] also lists constraints from combining the results from 
the WL probes with priors that CMB temperature and E-mode po- 
larization correlation measurements with an experiment like Planck 
could provide for the cosmological parameters other than /nl (see 
Appendix |A] for a computation of the prior Fisher matrix). How- 
ever, the improvement on the marginal uncertainties on /nl pre- 
dicted from this combination is very small. The reduction expected 
for the uncertainties in the other parameters is not large either ex- 
cept for the spectral index n^. 

To assess the robustness of our estimated parameter con- 
straints to changes of the realizations entering the covariance, we 



also compute the cosmological parameter constraints using the data 
covariance estimated from the 128 simulated fields of the fiducial 
zHORIZON runs. The resulting constraints are larger or smaller 
by up to few tens of percent, with a slight trend toward larger un- 
certainties for the results based on the zHORIZON covariances. 
For example, the marginal constraints from combining all three 
WL probes become: A/nl = 50, Afi,„ = 0.0018, Atrg = 0.0038, 
A/Is = 0.0084, Awo = 0.028. In comparison with the results listed 
in Table |2] the largest difference (a 25% increase) is observed for 

4.6 Biases induced by neglecting non-Gaussianity 

The constraints on the cosmological parameters predicted for a 
future all-sky weak lensing survey imply substantial correlations 
between /nl and other cosmological parameters. This raises the 
question whether neglecting a possible non-Gaussianity at the level 
of current constraints on /nl might induce a serious bias in the 
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Figure 12. Joint marginal constraints on pairs of cosmological parameters from a large future weak lensing survey: posterior means (symbols) and 68% 
confidence contours (lines) for the mean matter density O.^, the normalization erg and the spectral index fij, and the equation-of-state parameter wq obtained 
for a cosmology with 'true' non-Gaussianity parameter /ml = 30 when either allowing for /nl to vary (+, solid lines), or assuming /ml = (x, dashed lines) 
in the joint WL analysis. 



Table 3. Biases in cosmological parameters induced by assuming the 
'wrong' value for the primordial non-Gaussianity /nl in the analysis of a 
large future weak lensing survey: dilferences Afln,, Acrg, Aos, and Awo be- 
tween the expected posterior means and the 'true' values of the mean matter 
density the normalization erg, the spectral index n,, and the equation- 
of-state parameter »'o, resp., relative to the difference A/nl between the 
'true' value of the primordial non-Gaussianity /nl and its value assumed in 
the analysis. 
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other cosmological parameters. Such a bias would ruin the deter- 
mination of the cosmological parameters if the shift in the inferred 
posterior parameter confidence regions induced by neglecting non- 
Gaussianity were comparable to or larger than the inferred param- 
eter uncertainties. 

We study this possible bias on the basis of Eqs. l |32| l to l |34[ l, 
assuming flat priors for the cosmological parameters and the likeli- 
hood computed from our simulations. Using Eq. J33j, we can com- 



pute how the inferred posterior means of the other cosmological pa- 
rameters respond to a 'wrong' choice (i.e. assuming a value that dif- 
fers from the 'true' value) for /nl. This response is listed in Table|3] 
as ratio between the induced systematic error in the cosmological 
parameters and the systematic error A/nl in the non-Gaussianity 
parameter. 



For illustrative purposes, we consider a scenario where the 
'true' non-Gaussianity /nl = 30 (which is roughly the current 
best guess, [Komatsu et al.|2011^ . The joint marginal posterior con- 
straints on the parameters D,^, erg, n^, and vvo resulting from a joint 
bM are shown in Fig. 



12 



for two cases. In 



analysis of (A/^p), and < 
the first case, /nl is considered a parameter to be constrained from 
the data along with the other parameters. For our simple model, the 
resulting posterior means coincide with the 'true' values of the pa- 
rameters. The second case assumes /nl = in the analysis. This 
results in different means and covariances inferred for the parame- 
ters. The resulting biases are noticeable compared to the parameter 
uncertainties estimated for the full model including /nl. Moreover, 
the analysis ignoring a possible non-Gaussianity yields not only 
biased estimates, but also far too optimistic uncertainties for the re- 
maining parameters, in particular for erg and wq. Thus, the parame- 
ter biases are comparable to the inferred uncertainties in that case. 
As a result, the confidence regions inferred by an analysis ignoring 
non-Gaussianity do not contain the 'true' values of the cosmologi- 
cal parameters. 
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5 SUMMARY AND DISCUSSION 

Measurements of the level of non-Gaussianity in the primordial 
matter density fluctuations provide valuable information for con- 
straining models of inflation and the origin of cosmic structure. 
Currently the best constraints on primordial non-Gaussianity are 
from CMB experiments probing matter fluctuations at very high 
redshifts. However, lower-redshift observables could also provide 
competitive constraints. 

In this paper, we investigated how well a large weak lensing 
survey probing the matter distribution at redshifts < z < 1 , can 
constrain the level /nl of primordial non-Gaussianity of the local 
type. We employed a large set of A'-body simulations of structure 
formation in flat wCDM cosmologies, through which we performed 
full gravitational ray-tracing to simulate weak lensing observations 
similar to those proposed for Euclid and LSST. With these WL sim- 
ulations, we studied the influence of the cosmological parameters 
/nl, i^mi CTg. and wq on the cosmic shear correlation functions <f+ 
and the third-order aperture mass statistics {Ml^), and the shear 
peak abundance 0^^. Furthermore, we quantified the uncertainties 
and covariances of these WL probes, and the resulting constraints 
on the cosmological parameters. 

We found that if all other cosmological parameters were 
known accurately, each of the considered WL probes could pro- 
vide constraints A/nl ~ 20 on the non-Gaussianity parameter /ml. 
The constraints were much weaker (A/nl ~ 100 - 500), if the other 
cosmological parameters had to be constrained using the same WL 
data. However, the uncertainties in /nl could be reduced by com- 
bining information from the different WL probes. 

A joint analysis of the WL probes also yields much smaller 
uncertainties for the other cosmological parameters than may be 
obtained from any of the individual WL probes. When combining 
the probes, the expected marginal uncertainties at 68% confidence 
level on the cosmological parameters were: A/nl ~ 50, Aflni ~ 
0.002, Ao-g ~ 0.004, A/Is ~ 0.007, and Awo ~ 0.03. 

At this level of accuracy, neglecting non-vanishing non- 
Gaussianity could induce a noticeable bias in the determination of 
other cosmological parameters, in particular for the matter power 
spectrum normalization erg or the DE EOS parameter wq. More- 
over, neglecting a possible non-Gaussianity in the analysis may 
lead to far too optimistic estimates for the parameter errors. 

Our results for the expected conditional errors on /nl from 
cosmic shear are in agreement with the findings of |Fedeli and| 
|MoscardinT| j2010^ , who estimated errors of a few tens consider- 
ing the weak lensing power spectra. Our results for the expected 
constraints from shear peaks are somewhat less optimistic than the 
errors predicted by Matu ri et al.|f2011^ . This might be attributed 
to our more realistic predictions for the peak function, which we 
measured directly from simulated WL maps. 

The analysis in Section |4.4| showed that a large part of the 
parameter constraints stems from the lensing signals measured on 
small and intermediate angular or mass scales. This highlights the 
importance of improving the understanding of structure formation 
on these smaller scales, which is required for accurate model pre- 
dictions. Moreover, an improved understanding of the shape noise 
is needed for accurate predictions of the shear peak statistics at 
these scales. 

Our error analysis in Section [4~2] showed that for large future 
WL surveys, the statistical errors on the cosmic shear correlations 
and the aperture statistics are cosmic-variance dominated. This im- 
plies that for these probes, reducing the shape noise, e.g. by increas- 
ing the survey depth, will not reduce the statistical errors much. 



An obvious extension of this work is to take tomographic in- 
formation into account. Tomography might not significantly im- 
prove the errors on /nl arising from two-point statistics like the 
convergence power spectrum (see [Fedeli and Moscardini|[20T0l l, 
but it could improve the constraints from the peak abundance. It 
will also significantly reduce the errors on other parameters, such as 
the equation-of-state parameters of the Dark Energy. Such a study 
probably requires more simulations to obtain reliable covariances 
for the large data vectors arising in WL tomography. 

One should also investigate more the constraints one may ob- 
tain from combining the WL probes with other probes of large scale 
structure (see Gi annantonio et al.|2012| ). This will require simula- 
tions that include predictions for optical properties of galaxies or 
X-ray properties of clusters. Finally, a comprehensive future analy- 
sis should also consider other types of primordial non-Gaussianity 
and other cosmological parameters impacting structure formation. 
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APPENDIX A: PLANCK FISHER MATRIX 
Al Computing the CMB matrix 

In computing the Planck Fisher matrix we follow the methodol- 
ogy described in Eisenstein et al. ( 1999) and for the specific im- 
plementation we follow Ta kada and Jain ( 2009) . We thus assume 
that the CMB temperature and polarization spectra constrain 9 pa- 
rameters, and for our calculations we set their fiducial values to be 
those from the recent WMAP7 analysis ^Komatsu et al.|20TT| l. The 
fiducial parameters are: dark energy EOS parameters wo - -1.0 
and Wo = 0.0; the density parameter for dark energy fluE = 0.7274; 
the CDM and baryon density parameters scaled by the square of 
the dimensionless Hubble parameter w^. = Cl^h- = 0.1125 and 
Wb = il^h^ = 0.0226 {h = //o/[100kms-'Mpc-']); and the pri- 
mordial spectral index of scalar perturbations Ws = 0.963; the pri- 
mordial amplitude of scalar perturbations = 2.173 x 10 '; the 
running of the spectral index a = 0.0; and the optical depth to the 
last scattering surface t = 0.087. Hence we may write our vector 
of parameters: 



: (wo,Wa,ClDE,Wc,Wi„T,n„As,a)\ 



(Al) 



Whilst the CMB does constrain /nl through the temperature and 
polarization bispectra, here we are only interested in assessing the 
potential of weak lensing to place constraints on /nl. We therefore 
extend the vector p to the 10-dimensional vector: 



p ^ p = (wo,.. .,a,/NL)'- 



(A2) 



[T999|: 
F 



The CMB Fisher matrix can be written as jEisenstein et al.| 



ZZ¥^Cov-'[Q,,Q,]f- 



(A3) 



where {X, Y] e [TT, EE, TE), where Cf xr is the temperature power 
spectrum, ee is the E-mode polarization power spectrum, Cf xE 
is the temperature-E-mode polarization cross-power spectrum. We 
have been conservative and assumed that there will be no signif- 
icant information from the Cf bb> the B-mode polarization power 
spectrum. We compute all CMB spectra using CAMB ( [Lewis et al.| 
|2000[ l and use the additional module for time evolving dark energy 
models ^Hu and Sawick i 20071. We include information from all 
multipoles in the range (2 < £ < 1500). 

The covariance matrices for these observables are: 



CoV [C^ 77, C^ 7x] 
COV [Cf,TT, Cf,EE] 
CoV [Cf jTi C^je] 
CoV [Cf,EE, Cf,EE] 



1 



/sky2f+l 

1 2 2 

^^^Q.TE 

1 2 



.TTl 



/sky 2^+1 

1 2 



CfjE [Cf,TT + Nf rr] 



[Ce,EE + N( ,ee1 



(A4) 



(A5) 



(A6) 



(A7) 
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Cov[Cf,EE,QTE] 
COV [Cf,TE, C^je] 



1 



/sky 21 + 1 
1 



C[JE [Cf.EE + Aff,EE] 



+ 1 



/sky 2{ 

(Q,EE + ^f,TT)(Cf,EE + N(_ee)] 



(A8) 



(A9) 



In the above /^y is the fraction of sky that is surveyed and usable 
for science, and we take /j^y = 0.8. The terms A^f,TT ™d ^cee de- 
note the beam-noise in the temperature and polarization detectors, 
respectively. These can be expressed as: 



YT 



[wttM/4„„(^)J 



(AlO) 
(All) 



where wtt = [ATSseam] and w^e = [AE^Beam] • The beam win- 
dow function has the form: 



„m = expp(^+l)trL„]; ^Tb, 



(A 12) 



For the Planck experiment we assume that we have a single fre- 
quency band for science (143 GHz), and for this channel the fol- 
lowing parameters apply ( [The Planck Collaboration|2006^ : angular 
resolution of the beam SBeam = 7.1' [FWHM]; the beam intensity is 
At = 2.2(rcMB/lK)LuK], Ae = 4.2 (rcMB/lK) [pK]. We take the 
temperature of the CMB to be T = 2.726K. 



A2 Transforming from CMB to Weak Lensing parameters 

In the lensing simulations, we have considered how the weak lens- 
ing signal depends on the parameters, {wq, f2„,, «s> o"8. /nlI- Let us 
rewrite our original 10-D parameter set in terms of a new 10-D pa- 
rameters set. Therefore let us consider the transformation: 

P = (vi'o,W„,nDE,H'c, Wi„T,fIs,As,a,/NL)' 

1-^9 = {wo,w„,n^,h,fi„T,n„crg,a,ffn^)\ 



(A13) 



Six of the parameters are unchanged from the original set. The re- 
maining four are related to the original parameters in the following 
way: 



i^m - 1 ~ ^DE 

Wb 



o-g 



l-^DE 



1/2 



(A14) 
(A15) 

(A 16) 
(A17) 



where P(k\p) is the matter power spectrum, which depends on pa- 
rameters p„, and where the real space, spherical top-hat filter func- 
tion has the form Wi^iy) = 3[siny - >'cosy]/y', with y = kR and 
R = 8Mpc. 

We now need to construct the matrix A, i.e. the partial deriva- 
tives of the old parameters with respect to the new parameters. 
From Eqns l |A14[ l-l (A17^ we have q^, = Gaip), however in order 
to perform the partial derivatives we actually require the inverse 
of these relations, i.e. p„ = G^\q). In some cases these inverse 
relations may easily be determined, e.g. Eq. JaT4]i. However, in 
other cases no analytic inverse exists, e.g. Eq. ( |A17^ . A simple way 
around this problem is through recalling the following: 



(A 18) 



Table Al. Covariance Cpr of the non-Gaussianity parameter /nl, the mean 
matter density the normalization erg, the spectral index its, and the 
equation-of-state parameter wo assumed for a prior parameter distribution 
provided by a CMB experiment like Planck. 



/nl 



0"8 



"s 



/nl 

«s 



oo 

2.2x10-2 

-3.5x10-2 

-1.2x10-5 

1.0x10-' 



Or in other terms, 

'dp,. 





-3.5x10-2 
5.9x10-2 
-8.9x10-'' 
-1.7x10-' 



-St 





-1.2x10-5 l.OxlQ-' 

-8.9x10-* -1.7x10-' 

1.9x10-5 2.8x10-* 

2.8x10-* 5.8x10-' 



(A19) 



Aj^, we may then simply compute A, 



Hence, with dqa/dpy 
through: A = (A-')"'. 

Let us therefore form the matrix A-' 
untransformed parameters, we have A-j 
other four parameters, the situation is more complex and the non 
trivial elements are: 



For the case of the five 
: St- However, for the 



dSloE 
dA 
dwc 
dfi^ 
5wb 
dh 

di^DE 

dh 

dcrg 
dpa 



= -1 



Wb 



Wc + Wb 



[Wc -I- Wb]2 

[Wc + Wb]2 
1 

2 
1 
2 

1 

2o-g dpa 



{i-nuEy ' 

[(1 - nDE)(Wb + W + c)]-'/2 ; 

k dP{k\p) 



1 r d'/ 

2o-8 J (2n] 



f dpa 



WikRY 



(A20) 



(A21) 



(A22) 



(A23) 



(A24) 



(A25) 



Note, at leading order in /ml we have, dcrg/dfy^ = 0. In order 
to compute the derivatives dcrg/dpa we use the package CAMB. 
Besides the generation of various CMB power spectra, this pack- 
age can output the present day linear theory matter power spectra 
P{k\p). The derivatives are then determined numerically using the 
standard estimator for two sided derivatives. 

The numerical inverse of the matrix A-' can easily be com- 
puted using the SVD algorithm ( [Press et al.||1992| . We can then 
transform the CMB Fisher matrix from the original to the new pa- 
rameter set. The final step is a marginalization over those parame- 
ters not considered in our weak-lensing analysis. The covariances 
for the remaining parameters /nl, flm, o"8, "s, and wq are shown in 
Table IaT] 



APPENDIX B: PARAMETER COVARIANCE MATRICES 



Tables |B1|B8| show the posterior parameter covariances expected 
from the measurements of the cosmic shear correlation functions 
the aperture statistics {Mj^} and the shear peak abundance tpM 
in a large future weak lensing survey. Electronic versions of these 
tables are available online. 
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Table Bl. Posterior covariance of the non-Gaussianity parameter /nl, the 
mean matter density Cl^, the normalization erg, the spectral index Ms, and the 
equation-of-state parameter wq expected from measurements of the cosmic 
shear correlation functions f+. 



/nl 



£^8 



Vfo 



Al 

'h 
Wo 



4.5x10'* 

-5.1x10-' 
-2.5x10" 

-i.oxur' 

-l.OxlO' 



-5.1 xlO"' 

4.0x10-5 
3.0x10-5 
-1.4x10-"' 
6.9x10-" 



-2.5x10" 

3.0x10-5 
1.4x10-" 
-1.4x10-5 
6.2x10-" 



-l.OxlQ-' 
-1.4x10-" 
-1.4x10-5 

6.7x10-" 
-2.4x10"^ 



-l.OxlO' 
6.9x10-" 
6.2x10-" 
-2.4x10-3 
1.3x10-^ 



Table B2. Posterior parameter covariances expected from measurements 
of the third-order aperture statistics (M^p). 



/nl 



Wo 



/nl 

Wo 



2.0x105 
6.5x10-' 
-1.8x10' 
-l.lxlO' 

1.7x10' 



6.5x10-' 
1.0x10-" 
-1.7x10 



-4 



-4.0x10 



-4 



-1.8x10' 

-1.7x10-" 
1.8x10-3 
1.4x10-3 



1.6x10-3 -3.2x10-3 



-1.1x10' 

-4.0x10-" 
1.4x10-3 
2.1x10-3 

-6.9x10-3 



1.7x10' 

1.6x10-3 
-3.2x10 



-3 



-6.9x10 
3.5x10- 



-3 



Table B6. Posterior parameter covariances expected from joint measure- 
ments of the aperture statistics {M^^) and the shear peak abundance ^m- 



/nl 



0-8 



Wo 



/nl 

0-8 

«s 

Wo 



4.0x103 
1.8x10"^ 
-3.3x10-' 

-2.6x10-' 
-1.4x10" 



1.8x10-2 
1.6x10-5 
-2.0x10-5 
-4.5x10 



-5 



-3.3x10-1 
-2.0x10-5 
5.3x10-5 
6.1x10-5 



-2.6x10-' -1.4x10° 
-4.5x10-5 9.0x10-5 
6.1x10-5 -4.2x10-* 



2.5x10- 



-1.2x10- 



9.0x10-5 -4.2x10-* -1.2x10-5 1.5x10-3 



Table B7. Posterior parameter covariances expected from joint measure- 
ments of the cosmic shear correlation functions f±, the aperture statistics 
(Map), and the shear peak abundance ^m- 

/nl 

/nl 2.8xl03 



n„ 

0-8 
Wo 



-5.7x10-2 
-1.6x10-' 
1.1x10-' 
-1.2x10" 





0-8 




Wo 


5.7x10-2 


-1.6x10-' 


l.lxK)-' 


-1.2x10" 


2.8x10-* 


1.7x10-* 


-6.7x10-* 


3.5x10-5 


1.7x10-* 


1.3x10-5 


-7.3x10-* 


5.9x10-5 


■6.7x10-* 


-7.3x10-* 


4.5x10-5 


-5.7x10-5 


3.5x10-5 


5.9x10-5 


-5.7x10-5 


6.9x10-" 



Table B3. Posterior parameter covariances expected from measurements 
of the shear peak abundance <I>m- 





/nl 




0-8 


n 




Wo 


/nl 


4.8xl03 


l.OxK)-' 


-5.6x10-' 


2.3x10- 


-1 


3.5x10-' 




1.0x10-' 


7.0x10-5 


-1.1x10-" 


4.3xKr 


-5 


8.2x10-" 


0-8 


-5.6x10-' 


-1.1x10-" 


2.2x10-" 


-1.6xKr 


-4 


-1.4x10-3 


Is 


2.3x10-' 


4.3x10-5 


-1.6x10-" 


l.OxKr 


-3 


2.4x10-3 


Wo 


3.5x10-' 


8.2x10-" 


-1.4x10-3 


2.4x 10" 


-3 


1.4x10-2 



Table B4. Posterior parameter covariances expected from joint measure- 
ments of the cosmic shear correlation functions f ± and the aperture statistics 





/nl 




o-g 


"s 


Wo 


/nl 


2.0x10" 


-5.6x10-' 


-l.lxlO" 


1.4x10" 


-l.lxlO' 




-5.6x10-' 


1.8x10-5 


3.1x10-5 


-4.8x10-5 


3.2x10-" 


0-8 


-1.1x10" 


3.1x10-5 


6.8x10-5 


-8.5x10-5 


6.0x10-" 


Is 


1.4x10" 


-4.8x10-5 


-8.5x10-5 


1.6x10-" 


-8.2x10-" 


Wo 


-1.1x10' 


3.2x10-" 


6.0x10-" 


-8.2x10-" 


5.9x10-3 



Table B8. Posterior parameter covariances expected from joint measure- 
ments of the cosmic shear correlation functions f±, the aperture statistics 
(Map), and the shear peak abundance 0m> when combined with prior infor- 
mation from a CMB experiment like Planck. 



/nl 

Qui 
0-8 
"s 
Wo 



/nl 




0-8 




Wo 


2.6x103 


-4.5x10-2 


-1.4x10-' 


3.3x10-2 


-l.lxlO" 


-4.5x10-2 


2.1x10-* 


9.0x10-' 


-2.0x10-* 


2.9x10-5 


-1.4x10-' 


9.0x10-'' 


1.2x10-5 


-2.2x10-* 


5.1x10-5 


3.3x10-2 


-2.0x10-* 


-2.2x10-* 


1.3x10-5 


-1.7x10-5 


-1.1x10" 


2.9x10-5 


5.1x10-5 


-1.7x10-5 


6.4x10-" 



Table B5. Posterior parameter covariances expected from joint measure- 
ments of the cosmic shear correlation functions f ± and the shear peak abun- 
dance ^M- 





/nl 




0-8 




wo 


/nl 


3.6x103 


-8.2x10-2 


-2.0x10-' 


2.2x10-' 


-1.4x10" 


Qni 


-8.2x10-2 


3.9x10-* 


2.5x10-* 


-1.1x10-5 


4.5x10-5 


0-8 


-2.0x10-' 


2.5x10-* 


1.5x10-5 


-1.2x10-5 


6.3x10-5 


Is 


2.2x10-' 


-1.1x10-5 


-1.2x10-5 


7.1x10-5 


-8.6x10-5 


Wo 


-1.4x10" 


4.5x10-5 


6.3x10-5 


-8.6x10-5 


7.7x10-" 
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